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ABSTRACT 

A  pair  of  complex  conjugate  fourth-order  differential  equations  that 
govern  the  deformation  of  orthotropic  circular  cylindrical  shells  is 
presented.  As  shown  in  the  paper,  this  pair  of  equations  is  as  accurate  as 
equations  can  be  within  the  scope  of  the  Kirchhoff  assumptions.  Also 
presented  for  the  first  time  are  several  pairs  of  accurate  and  simple  fourth 
order  equations  which  can  be  systematically  and  explicitly  deduced  from  the 
previously  mentioned  pair  of  equations.  Because  of  their  accuracy  and 
simplicity,  these  simple  equations  are  of  practical  importance.  The  advantage 
in  applying  those  equations  presented  herein  is  that  their  solutions  can  be 
easily  found  in  simple  closed  forms.  This  considerably  simplifies  calcula¬ 
tions  for  solving  problems  of  orthotropic  and  laminated  composite  cylindrical 
shells.  Unlike  other  known  equations  in  the  literature,  their  general  solu¬ 
tions  remain  unknown  because  of  the  algebraic  complexities  involved. 
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AN  ACCURATE  THEORY  AND  SIMPLE  FOURTH  ORDER  GOVERNING 
EQUATIONS  FOR  ORTHOTROPIC  ANC  COMPOSITE  CYLINDRICAL  SHELLS 

Shun  Cheng*  and  F.  B.  He* 

Introduction 

Considerable  attention  has  been  devoted  to  the  study  of  Isotropic  cylin¬ 
drical  shells.  Literature  on  thts  subject  is  quite  extensive.  In  contrast, 
relatively  little  work  has  been  done  on  the  formulation  of  the  basic  equa¬ 
tions  for  orthotropic  cyl Indr teal  shells,  although  they  are  frequently 
employed  as  structural  elements  In  Industry  [1-7].  Examples  of  orthotropic 
shells  Include  laminated  composite,  perforated  and  stiffened  cylindrical 
shells  whose  material  behavior  can  be  considered  as  orthotropic.  Composite 
shells  (5,6,7)  constitute  an  example  of  great  practical  Importance. 

As  is  known  in  the  literature,  the  classical  shell  theory  is  based  on 
the  same  basic  assumptions  emplo^d'  in  the  theory  of  thin  plates,  known  as 
Klrchhoff  assumptions.  Since  the  inception  of  Love's  first  approximation, 
further  simplifications  or  approximations  beyond  these  basic  assumptions 
have  been  Introduced  In  developing  the  theory  of  thin  shells.  As  the  abundance 
of  literature  Indicates,  many  versions  of  shell  theories  have  been  formulated, 
each  depending  on  different  versions  of  the  various  approximations.  This  has 
confronted  engineers  as  well  as  researchers  with  a  controversial  problem  with 
regard  to  the  consistency  of  the  theory,  shortcomings  of  the  derivations  and 
accuracy  of  the  resulting  equations.  Many  sets  of  resulting  equations  have 
been  proposed  for  Isotropic  shells  [8-12]  and  especially,  due  to  their 
Importance  In  application  and  the  fact  that  they  display  nearly  every  type 
of  behavior  found  In  general  shell  theory,  for  cylindrical  shells  [8-16]. 
Several  publications  are  listed  In  the  reference.  Others  may  be 
found  by  consulting  these  references.  As  for  orthotropic  cylindrical  shells, 
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two  types  of  basic  equations,  corresponding  either  to  Flugges's  or  Donnell's 

equations  for  Isotropic  shells,  have  been  formulated  in  the  literature  [1-3,6]. 
In  either  case,  a  resulting  single  eighth-order  differential  equation  maybe 

deduced.  However,  the  eighth-order  equation  for  orthotropic  shells  is  more 
complicated  than  the  corresponding  ones  for  isotropic  shells.  A  common 
difficulty  with  these  eighth-order  equations  in  isotropic  or  orthotropic 
shell  theory  is  that  their  general  solutions  remain  unknown  because  of  the 
algebraic  complexities  Involved.  For  orthotropic  cylindrical  shells,  even 
the  simpler  eighth-order  equation  based  on  Donnell  approximations  as  seen  in 
[2]  suffers  from  the  same  complexity.  Although  a  fourth-order  equation  is 
presented  In  [3],  this  equation  does  not  yield  accurate  and  dependable 
solutions  as  is  Illustrated  through  the  comparison  made  between  numerical 
results  from  analytical  solutions  and  from  experimental  data  presented  in 
the  same  paper  [3].  In  computing  the  characteristic  roots  arising  from 
solving  these  eighth-order  equations  by  means  of  eigenfunctions,  it  is  found 
that  the  two  large  roots  and  the  two  small  roots  in  the  same  set  of  solutions 
for  the  characteristic  equation  are  far  apart  and  of  different  orders  of 
magnitude.  This  makes  the  computation  more  tedious  and  time-consuming, 

‘even  with  the  present  day  numerical  techniques. 

Recently  a  general  theory  for  thin  Isotropic  shells  was  developed  by 
Markov  [17].  It  Is  a  consistent  theory,  since  it  makes  no  simplifications 
or  approximations  beyond  a  clear  set  of  fundamental  hypotheses.  Other 
advantages  of  the  method  of  derivation  as  applied  to  shells  of  general 
curvature  have  also  been  illustrated  in  [17]. 

In  the  present  paper,  a  pair  of  complex  conjugate  fourth-order  partial 
differential  equations  that  govern  the  deformation  of  orthotropic  circular 
cylindrical  shells  is  proposed.  This  pair  of  equations  is  deduced  from  a 


set  of  basic  equations  which  Is  based  on  the  following  Kirchhoff  hypotheses: 
(a)  The  transverse  normal  stress  is  negligibly  small  and 


(b)  Normals  to  the  middle  surface  of  the  shell  remain  normal  to  it  and 


undergo  no  change  in  length  during  deformation. 

The  set  of  basic  equations  is  exact  in  the  sense  that  in  deriving  these  equa¬ 
tions  all  terms  have  been  retained  without  introducing  further  simplifications 
or  approximations  beyond  these  fundamental  hypotheses.  Even  those  terms  which 
are  of  higher-order  are  kept  since  they  can  be  summed  in  closed  form. 

Because  the  pair  of  equations  deduced  herein  is  complex  conjugates, 
only  one  of  the  equations  needs  to  be  considered.  Further,  closed  form  solu¬ 
tions  of  the  characteristic  equations  that  arise  from  solving  the  pair  of 
governing  equations  by  means  of  eigenfunctions  can  be  easily  obtained.  The 
technique  used  is  an  extension  of  the  one  for  isotropic  shells  presented  in 
[15,16].  From  the  pair  of  equations,  a  number  of  simplified  fourth-order 
governing  equations  can  be  systematically  and  explicitly  deduced,  as  shown 
in  the  paper.  These  fourth-order  equations  for  orthotropic  cylindrical  shells 
are  new  in  the  literature  and  of  definite  technical  importance  because  these 
equations  can  be  easily  solved  in  closed  forms  and  yet  retain  practically 
the  same  accuracy  as  the  original  eighth-order  equation. 

Basic  Equations 

In  accordance  with  the  fundamental  hypotheses  stated  previously,  the 
following  basic  equations  can  be  deduced  for  orthotropic  circular  cylindrical 
shells.  Let  a  be  the  radius  of  the  midsurface  of  the  shell,  x,  y,  z  the 
axial,  circumferential  and  radial  coordinates  and  a,  3  the  dimensionless 
midsurface  coordinates  along  lines  of  curvatures  (a  ■  -  ,  3  *  ^).  The  three 

d  d 

displacement  components  u^,  u^,  and  u2  of  an  arbitrary  point  of  the  shell 
can  be  expressed  in  terms  of  midsurfae  displacements  u,  v,  and  w  as 
follows  [8,16] 


The  stress-strain  relations  for  orthotropic  materials  [18,1 9]  are 

E1  E2 
°o  '  <  We*  •  °$  ■  <  We1  1 
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where  Ej ,  E2  are  the  moduli  of  elasticity  along  the  principal  directions 
a  and  8,  respectively,  G  Is  the  shear  modulus  which  characterizes  the 
change  of  angles  between  principal  directions  a  and  0,  Is  the 

Poisson's  ratio  which  characterizes  the  decrease  in  a-direction  due  to 
tension  applied  In  6-dlrectlon,  and  v2  *  is  the  Poisson's  ratio  which 
characterizes  the  decrease  In  8-dlrectlon  due  to  tension  applied  in  a-direc 
tlon.  Among  these  material  constants  there  exists  the  relation  [18,19]: 

Elvl  *  E2v2  (3) 

‘The  components  of  strain  at  an  arbitrary  point  of  the  shell  are  related 
to  the  midsurface  displacements  by  [8,15,16] 
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The  bending  and  twisting  (t)  strains  are 
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Let  h  be  the  wall  thickness,  K1 ,  K2  the  extensional  rigidity,  ,  02  tl 
flexural  rigidity 
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Let  N  ,  N-  be  the  normal  stress  resultants,  S  .  Sfl  the  shear  stress 
resultants.  Mo»  the  bending  moments,  Mag,  the  twisting  moments,  and 
Qp  the  transverse  stress  resultants  [15].  These  are  stress  resultants 
(N,  S,  Ql  and  couples  (M)  per  unit  length  of  the  middle  surface  and  are 
related  to  the  midsurface  displacements  through  the  stress-strain  relations 
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In  which 
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The  equations  of  static  equilibrium  are 

3N  as  3N  as 
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In  which  X,  Y  and  Z  are  surface  loads  per  unit  area  in  x,  y  and  z 
directions,  respectively. 


Pair  of  Accurate  Complex  Conjugate  Fourth-Order  Equations  for  Normal  Deflection 
Substituting  equation  (8)  into  equation  (9),  a  system  of  three  differ¬ 
ential  equations  is  obtained  for  the  three  basic  functions.  This  system  is 
presented  In  Table  1  and  possesses  a  symmetrical  structure.  The  three  linear 
partial  differential  equations  with  constant  coefficients  can  be  reduced  to 
a  single  differential  equation  of  higher  order  that  is  more  convenient  to 
solve  and/or  analyze  with  regard  to  the  present  problem.  These  three  equations 


presented  In  Table  1  shall  be  considered  as  algebraic  equations  in  u,  v,  and 

o 

w  having  coefficients  which  are  constants  (elastic  constants  and  c  )  and 
the  symbols  of  differentiations.  Let  DQ  be  the  3x3  determinant  of  Table  1 
and  calculate  Its  cofactors  D^,  D^.  ....  D33<  Let 

u  *  D11*1  *  v  a  °i2<*>1  •  w  *  d13^i  (sum  on  1  •  1  a  1 .2,3)  (10) 

and  substitute  these  expressions  in  the  three  equations  in  Table  1.  Then,  in 
accordance  with  the  theory  of  linear  algebra 

1  "  vlv2  2 

Vi  *  "  e/"  a  X1  *  0  »  0  *  1.2.3)  (11) 

are  obtained,  in  which  X1  =  X,  Xg  ®  T,  X3  *  Z.  If  only  a  normal  surface 
load  Z  is  applied  on  the  shell,  ^  and  <f>2  can  be  set  equal  to  zero  in 

equations  (10)  and  (11).  Calculating  cofactors  D31 ,  D^,  D33  and  Dq  from 

Table  1  and  replacing  <J>3  by  ■£-  ,  the  following  are  obtained  from 
equations  (10)  and  (11) 
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The  constant  coefficients  In  equations  (121,  03)  and  04)  may  contain  co- 

p 

efficients  higher  than  those  shown  (of  the  order  c  or  higher).  These 

coefficients  have  been  omitted  since,  tn  thin  shell  theory,  j  <  jjy  »  thus 
2  -4 

c  <  2x10  Is  a  very  small  number.  The  complete  expression  of  DQ  Is 
given  in  the  Appendix.  Comparison  of  the  magnitude  of  the  coefficients 
of  the  terms  which  were  omitted  with  the  coefficients  of  the  terms  (having 
the  same  partial  differentiations)  which  were  retained  in  equation  (14) 
reveals  that  these  omitted  terms  are  truly  of  smaller  orders  of  magnitude. 
This  fact  has  been  further  verified  through  the  actual  computation  of  thes 
coefficients  using  available  numerical  data  drawn  from  elastic  constants  o 
many  orthotropic  materials.  In  Table  2,  the  elastic  constants  of  a  few 
materials  are  presented.  Thu s# equation  (14)  is  an  accurate  governing 
equation  for  orthotropic  cylindrical  shells  because  this  equation  is  derived 


from  the  basic  hypotheses  without  introducing  further  approximations  in  its 
derivation  except  that,  as  just  stated,  some  negligibly  small  terms  have  been 
omitted.  These  small  terms  have  been  totally  dropped  in  all  the  known  equa¬ 
tions  of  orthotropic  shells.  In  some  publications,  even  certain  terms  in 
equation  (14)  are  neglected.  In  the  following  analysis,  some  of  these 
negligibly  small  terms  will  be  retained  so  that  equation  (14)  can  be  reduced 
to  a  pair  of  fourth  order  complex  conjugate  equations.  This  not  only 
tremendously  simplifies  calculation  of  the  roots  of  the  characteristic 
equation  which  arise  from  solving  the  equation  by  separation  of  variables 
but  It  also  facilitates  obtaining  solutions  in  simple  explicit  forms.  As 
stated  previously,  finding  solutions  of  equation  (14)  in  explicit  forms  is 
almost  prohibitively  difficult  due  to  the  algebraic  complexity  involved. 

In  addition  to  keeping  some  small  terms,  the  following  approximate  relation  as 


given  In  [2Q,21]  Is  employed: 
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2(1 
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Table  2  Mechanical  Properties  of  Materials 


Species 

k  =  E2/E1 

G/Ej 

V1 

v2 

Glass/epoxy 

0.3333 

0.1666 

0.0833 

0.2500 

Boron/epoxy 

0.1000 

0.0333 

0.0300 

0.3000 

Graphite/epoxy 

0.0250 

0.0125 

0.0063 

0.2500 

Douglas-fir 

0.0500 

0.0780 

0.0220 

0.4490 

Table  3  Characteristic  Roots  (w 
(Boron/Epoxy) 


cosna) 


Equation 


50  0.001 


50  0.010 


0.009297 

0.009278 

0.009278 

0.008900 

0.008872 

0.008872 

0.103515 

0.103496 


0.093577 

0.093411 

0.093414 

0.089696 

0.089463 

0.089420 

0.328262 

0.327653 


+  Q.008429i 
+  0 . 00841 4i 
+  0.00841 4i 
+  0. 008843 i 
+  0.008871 i 
+  0.008871 i 
+  0.042850i 
+  0. 042858i 


+  0. 083514i 
+  0. 083407 i 
+  0.083402i 
+  0.0877361* 
+  0.088064i 
+  0. 088024 i 
+  0.1 351 23i 
+  0 . 1 35377 i 


0.100 


3.217917 

3.217926 

3.211572 

3.217595 

3.211551 

3.282622 

3.276528 


+  1 . 37697 Oi 
+  1.376965i 
+  1 .3801481 
+  1.3777901 
+  1 .380201  i 
+  1  .3512341 
+  1.3537691 


P3.  P4 


0.00°079 

0.000079 

0.000079 

0.000079 

0.000079 

0.000079 

0.042877 

0.042869 


0.007852 

0.007857 

0.007857 

0.007870 

0.007883 

0.007872 

0.135969 

0.135719 


1.331 
1.331469 
1.331472 
1.329567 
1 .332245 
1 .329618 
1.359694 
1.357185 


+  1 .0000001 
+  1 .0000001 
+  1 .OOOOOOi 
+  0. 999997 i 
+  0.999992i 
+  0. 999997 i 
+  0.1 03450i 
+  0.1034691 


+  1 . 0001 391 
+  1  . 0001 38 1 
+  1 . 0001 38i 
+  0.999889i 
+  0. 9993721* 
+  0.999839i 
+  0.3262201 
+  0.326828i 


+  3 .327 954 i 
+  3 .327  934i 
+  3 .327  925i 
*  3 .3337501 
+  3 . 327594 i 
+  3.3337281 
+  3.2621 971 
+  3 .2682801 
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11. 


in  which  and  ,^V*2  are  9eometr*c  mean  values  for  the  modulus  E 
and  Poisson's  ratio  v,  respectively.  From  equations  (3),  (15)  and  (16), 
we  obtain 


2kj  »  vTTO-v^vJ),  K  -  2(2k]+v1)=2/IT,  k£  =  6k  (17) 

a8 

In  equation  (14),  replacing  the  coefficient  k9  of  the  middle  term  — j — r 

c  3a  33 

(only  in  this  term)  by  6k  and  keeping  some  of  these  omitted  small  terms, 
from  equations  (13)  and  (14),  the  governing  differential  equation  for  normal 
deflection  w  may  be  written  as 


_  .4  -4  .  ?  jv4  «4 

L  L  w  =  ^  (k-2v.krv?)  -§-*■  +  k  -ix]  Z 

U1  3aS  k1  1,1  3a  38  33 


(18) 


in  which 


4  4  4  2  7 

.  .  3  +  „  3  *  wf  3  *  3  *  3  ^ 

L - r  +  K  — a — t  +  J  +  — T  + 

3a  3a  dp  33  3a  33 

1  5^  ^4  ,.2  ~2 

+  i  [- - o"  +  c-i(k~  — j-  +  k*  — s — k  +  kc  — j"  +  — 9  +  kfi  — 5)] 

C1  3a^  1  J  3a*  *  3a^33  33  5  33^  5  3a* 

(19) 


and  L  is  the  complex  conjugate  linear  differential  operator  of  L,  i  =  and 


C1  =  . ;  »  k.  =  v,  -  k  ,  k4  =  1  (k2-2kK-2k)  , 

/k(l-v1  \^2 ) 

k5  =  k(K-k-v^ )  ,  kg  =  (4k-k2-3v2) 


(20) 


As  will  be  shown  later,  the  replacement  of  coefficient  k2  of  the  middle 
term  in  equation  (14)  by  6k  has  only  a  negligibly  small  effect  on  the  solu¬ 
tions  of  the  problem,  although  the  expression  (16)  is  merely  an  approximate 
relation  and  may  not  be  as  accurate  as  is  the  corresponding  relation  for 

isotropic  materials.  In  equation  (19),  it  is  seen  that  the  last  term 

a2  2 
C1  0  *s  very  small  as  compared  with  the  term  —  -2-=- 

3<>  C1  3a2 

hence  can  be  dropped.  This  should  not  yield  any  noticeable 


and 


effect  on  the  accuracy  of  the  equation  as  will  be  further  elaborated  later. 

The  homogeneous  solutions  of  equation  (18)  are  obtained  from 

Lw  =  0,  Cw  »  0  (21 ) 

From  equations  (12)  and  (13),  we  can  express  u  and  v  in  terms  of  w  [15]. 
Equation  (18)  reduces  to  the  same  governing  equation  for  isotropic  cylindrical 
shells  as  deduced  in  D9.  Using  the  relation  (16),  equation  (18)  can  also  be 
reduced  to  the  orthotropic  plate  equation  [19]  as  the  radius  of  the  shell 
goes  to  infinity. 

Solutions  by  Eigenfunctions 

It  may  be  shown  that  homogeneous  equation  (21)  and  suitable  boundary 
conditions  are  satisfied  by  making  use  of  the  following  solution  when  the 
eigenfunctions  are  trigonometric  along  a  generator: 

w  *  ep6cosna  (22) 

in  which  n  =  ,  m  is  an  arbitrary  integer,  £  represents  the  length  of 

the  shell  and  e  is  the  base  of  natural  logarithms.  When  the  eigenfunctions  are 
trigonometric  in  the  circumferential  direction,  w  can  be  taken  as 

w  =  epctcosn8  (23) 

in  which  n  is  a  real  number.  It  is  an  integer  value  when  the  cylinder  is 
closed  and  a  noninteger  value  when  the  shell  is  open.  Substituting  expres¬ 
sions  (22)  and  (23)  into  the  governing  equation  (21)  yields  characteristic 
equations  for  the  determination  of  the  roots  p.  Four  complex  roots  are 

obtained  and  the  other  four  roots  are  the  complex  conjugate  numbers  to  these 

2 

four  roots.  The  characteristic  equations  are  quadratic  equations  in  p  . 


Hence  solutions  of  the  present  problem  can  be  easily  found  in  closed  forms. 


Simple  Equations 

The  accurate  fourth-order  equation  (21)  can  be  used  to  obtain  a  number 

of  simplified  equations  which  are  new  In  the  literature  and  are  of  Importance 

In  practice.  Considering  the  actual  values  of  elastic  constants  of  various 

2 

orthotropic  materials  and  the  smallness  of  c  ,  It  can  be  easily  shown  that 

a2 

the  last  term  in  equation  (19)  c-jkg-^-j  Is  much  smaller  than  the  term 
l  a2 

- — ,  hence  this  term  can  be  dropped  In  equation  (19)  as  previously  stated. 
C1  3a 

When  the  same  considerations  are  applied,  terms  with  coefficients  k^,  k^  and 
kg  In  equation  (19)  can  also  be  neglected  because  they  are  of  a  smaller  order 
of  magnitude  In  comparison  with  other  terms  which  have  the  same  partial 
differentiations  In  the  equation.  Dropping  these  terms  in  equation  (19) 
yields  the  following  simplified  equation 


Lw  =  [• 


+  K 


da 


34  ,  k  34  ,  |c/  32  ,  32 

3?  V  V 


)  ±  J-  -4]  w  =  0 

C1  3a^ 


(24) 


If,  in  equation  (24),  new  dimensionless  coordinates  £  and  c  are  introduced 

by  stretching  the  variables  a  and  0  such  that  a=/c,  £  and  6  =  n. 

1  2  2 
a  d 

the  fourth  term  and  the  fifth  term  in  equation  (24)  become  kc, ( — s-  +  — 7). 

.  Z  dC  3n 

The  fourth  term  is  small  as  compared  with  the  last  term  - — and  hence 

C1  3C 


has  little  effect  on  the  characteristic  roots.  Therefore  this  term  can  be 
dropped  in  equation  (24)  and  another  simple  equation 


Lw  ■  (- 


+  K 


3a 


2  1 

3a4  3  3 


+  k 


38  33  1  3ct 


(25) 


is  obtained.  If  the  fourth  term  in  equation  (25)  is  also  dropped,  one 
obtains 


Lw 


.4* 


— 7 — 7 


+  k 


a4  i  a2 
3  *07^"' 


(26) 


14. 


Unlike  others  presented  in  this  section,  the  preceding  equation  is 
obtained  by  dropping  a  term  in  equation  (25)  without  the  usual  presence  in 
that  equation  of  a  similar  term  having  a  much  larger  coefficient  with  which, 
as  a  justification  for  the  resulting  simplification,  the  dropped  term  can  be 
compared.  This  procedure  unavoidably  affects  the  accuracy  of  the  equation, 
as  will  be  seen  in  the  next  section. 

In  deducing  the  fourth-order  equation  (21)  from  the  original  equation  (14) 

a8 

only  the  coefficient  k9  of  the  term  — t— - *  in  equation  (14)  is  replaced  by 

6  3a  36 

its  approximate  value  6k  (17).  If  the  same  approximation  is  also  used  for  the 

coefficient  K  in  equation  (19),  then  equations  (24),  (25)  and  (26)  become 

2 

(V4  +  kv2  ±  -L  8  )  w  =  0  (27) 

0  C1  3ct 

(V*  +  k  -4  ±  r  -4)  w  -  0  (28) 

0  38  ci  3ci 


(ij  *£77)  ws0 


in  which 


v?  =  4* 


_2  ^2 
r  =  -S—  +  *  ,  k  =  ^  and  c, 

3a  36^  1 


”vlv2 


By  substituting  the  differential  operators  L  given  by  equations  (24),  (25) 
and  (26)  into  the  left-side  of  equation  (18),  the  complete  version  of  these 
equations  including  the  load  term  can  be  written  as 

L  c’"  1 2<K-2k’‘v< 1  Ji?  * k  i?3  2  (30) 


in  which  L  represents  any  one  of  the  three  linear  differential  operators  of 
equations  (24),  (25)  and  (26).  Similarly  the  complete  version  of  equations 
(27),  (28)  and  (29)  are 

L  Lw  =  jr-  V*  Z 

D-j  0 


(31) 


15. 


In  which  L  represents  any  one  of  the  three  linear  differential  operators 
of  equations  (27),  (28)  and  (29l).  The  Morley,  Novozhilov  and  Donnell  equa¬ 
tions  for  Isotropic  shells  [15]  are  special  cases  of  the  equations  (27),  (28) 
and  (29). 


which  can  be  further  simplified  to 


[ — j  +  i 
da4 


rT-vTvTT  h2 

1  Zi  r  d 

” ][7T  "  1 
da 


^(l-v1v2)  (1-v.  v2)h 

- - - ]w  =  - r  z 

C  12  E-jC^ 


Problems  of  Thin  Rings  and  Long  Tubes 

When  a  ring  is  loaded  by  forces  applied  at  the  boundary,  parallel  to 

the  plane  of  the  ring,  the  stress  components  are  zero  on  both  faces  of  the 

ring.  Such  a  state  of  stress  is  called  plane  stress.  Following  the 

•procedures  presented  in  [22] 

2  2 

7?  h' 

<34> 

where  I  *  bh^/12,  b  *  width  of  the  ring,  M  is  the  bending  moment  (M  *  bM2), 
and  u>  represents  the  rotation  of  radial  cross  sections  of  the  ring.  For 
an  infinitely  large  radius  a  the  preceding  equations  coincide  with  that 
for  a  straight  beam. 


When  a  long  circular  tube  is  under  the  action  of  lateral  loads  uniformly 
distributed  along  the  axis  of  the  cylinder,  we  have  a  state  of  plane  strain. 


In  this  case,  displacement  along  the  axis  of  the  tube  u  Is  zero  and  v 
and  w  are  functions  of  8  only.  Following  the  procedures  similar  to  the 
deductions  of  the  basic  equations  for  thin  rings,  the  following  basic  equa¬ 


tions  for  the  bending  of  long  tubes  of  orthotropic  materials  can  be  obtained: 

dV„  >2(1-V2la2M2 

dr  e0  tr 


dv 

38 


+  w  *  0 


(35) 


Thus  the  basic  equations  of  the  present  theory  contain  both  ring  bending 
and  bending  of  long  circular  tubes  as  special  cases.  However,  as  stated  in 
[16],  the  equations  for  bending  of  thin  rings  and  long  tubes  cannot  be 
deduced  from  the  Donnell  equations. 

Comparisons  and  Conclusions 

Utilizing  the  computer,  the  relative  accuracy  of  the  differential  equa¬ 
tions  presented  previously  can  be  further  studied  through  numerical  techniques. 
This  can  be  done  by  calculating  the  numerical  values  of  the  characteristic 
roots  of  the  equations  and  making  a  comparison  of  the  closeness  of  these 
roots.  The  elastic  constants  of  several  typical  orthotropic  materials  are 
collected  in  Table  2.  Using  these  values,  the  roots  calculated  from  homogeneous 
equations  of  (14)  and  (18)  and  equations  (24-29)  are  obtained.  Many  roots 
for  other  orthotropic  materials  have  also  been  calculated.  Similarities  in 
the  properties  of  these  roots  for  different  materials  can  be  observed.  However, 
due  to  space  limitations  and  the  fact  that  the  same  conclusions  can  be 

drawn  from  different  materials,  only  Boron-epoxy,  Glass-epoxy  and  Graphite- 
epoxy  are  presented  in  Table  3  through  Table  8  for  a  range  of  significant 
parameters.  From  all  the  numerical  results. 


17. 


It  may  be  concluded  that  the  differences  between  the  roots  of  equations  (14), 

(18)  and  (24)  are  negligibly  small  for  all  values  of  n  and  c.  Numerical 
results  also  show  that  simplified  equations  (25),  (27)  and  (28)  can  yield 
accurate  solutions  as  seen  from  the  closeness  of  the  characteristic  roots 
of  these  equations  to  those  of  equations  (14),  (18)  or  (24).  As  is  expected, 
the  thinner  the  shell,  the  closer  In  value  these  roots  will  be.  The  simpl1> 
fled  equations  (26)  or  (29),  which  Is  only  one  term  less  than  equation  (25) 
or  (28),  Is  not  always  as  accurate  and  dependable  as  other  equations  £16].  These 
two  equations  are  apparently  inaccurate  In  the  case  w  »  ep^cosna  when  n 

Is  small.  Hence,  special  care  is  needed  when  they  are  employed.  All  the  preceding 
conclusions  hold  also  for  the  case  when  E1  and  E^,  v1  and  are 

Interchanged  In  the  calculations  and  can  be  applied  to  laminated  composite 
shells.  In  conclusion,  equations  (24),  (25),  (27)  and  (28)  deduced  herein 
have  the  two  essential  properties  of  accuracy  and  simplicity  and  hence  are 
of  importance  in  appl Ications. 


.I 
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Table  4  Characteristic  Roots  (w  *  epacosn8) 
(Boron /Epoxy) 


1/C 

n 

pl»  p2 

P3.  P4 

(14) 

2.802685  +  2.8080321 

0 

(18) 

2.803801  +  2.8091361 

0 

(24) 

2.803801  +  2.8091361 

0 

50 

2.796462  +  2.8142841 

0 

(25) 

2.805359  +  2.8053591 

0 

(28) 

2.805359  +  2.8053591 

0 

(26) 

2.805359  +  2.8053591 

0 

2.805359  +  2.8053591 

0 

TO 

HV'JL’kL'IiT.'fin'JrnK.'VL'nB 

0 

(18) 

2.945199  +  2.6717301 

0 

(24) 

2.945332  +  2.6716091 

0 

50 

l 

(27) 

2.853203  +  2.7583171 

0 

(25) 

2.948208  +  2.6694311 

0 

2.862275  +  2.7495751 

0 

(26) 

2.948680  +  2.6700701 

0.058927  ♦  0.0533591 

(29) 

2.862805  +  2.7501731 

0.057446  +  0.0551861 

|B|HB 

(14) 

ag.WO  t  rt*n.W»)  UMI 

n 

(18) 

28.109206  +  27.9980891 

0.019553  +  0.0194951 

(24) 

28.109254  +  27. 998038 i 

0.019553  +  0.0194951 

0 

(27) 

28.075257  +  28.0319501 

0.019539  +  0.0195091 

C 

(25) 

28.109475  ♦  27.9978271 

0.019563  +  0.0194851 

9  K 

(28) 

28.076149  +  28.0310601 

0.019540  +  0.0195081 

(26) 

28.109477  +  27.9978291 

0.022589  +  0.0225001 

iHH 

_ m 

28.076151  +  28.0310621 

0.022563  +  0.0225261 
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Table  S 


Characteristic  Roots  {w  ■  ep®cosnot) 
(Glass/Epoxy) 


Equation 


0.006679 

0.006684 

0.006688 

0.0Q6585 

0.006546 

0.006546 

0.088900 

0.088900 


0.067138 

0.067147 

0.067149 

0.066126 

0.065748 

0.065748 

0.281398 

0.281363 


2.738576 

2.738575 

2.738617 

2.738139 

2.738452 

2.738106 

2.813983 

2.813632 


♦  0.0063 93 i 
4  0.0064 06 i 
+  0. 006403 i 

♦  0. 006508 i 

♦  0.0065461 

♦  0.0065461 

♦  0.03681 6i 

♦  0.0368171 


37291 

♦  0.063739i 

♦  0.0637411 

♦  0.0648011 
+  0.0651 87 i 

♦  0.0651871 

♦  0.11 6309i 

♦  0.1163271 


4  1.193  2  5  0i 
+  1.1932511 

♦  1.1932671 
4  1.1936491 

♦  1.1935651 
4  1.1937331 
4  1 .1 630891 
4  1.1632681 


p3*  P4 


0.000043 

0.000043 

0.000043 

0.000043 

0.000043 

0.000043 

0.036824 

0.036824 


0.004284 

0.004285 

0.004285 

0.004287 

0.004286 

0.116559 

0.116545 


1.132930 

1.132931 
1.132948 
1.132994 
1.133228 
1.133074 
1.165591 
1.165446 


4  1.0000001 
4  1 . 000000 i 
4  1 .0000001 

4  0.9999991 
4  0.9999981 
4  0. 9999981 

4  0.0888811 

4  0.0888331 


.0 

4  1.0000331 
4  1.0000331 
4  0.9999231 
4  0.99984 7 i 
4  0.9998731 
4  0.2807941 
4  0. 280837 i 


4  2.8843871 
4  2.8843871 
4  2. 8844*51 
4  2.8847261 
4  2.8842561 
4  2.8846911 
4  2.8079421 
4  2.8083731 


Table  6  Characteristic  Roots  (w  =  epacosn6) 
(Glass/Epoxy) 


5000  0.700 


0.013357  +  0.0117291 
0.013303  +  0.011 695i 
0.013304  +  0.0116941 
0.012591  +  0.0125481 
0.012572  +  0.0125681 
0.012571  +  0.0125671 
0.123248  +  0.0509921 
0.123199  +  0.0510121 

0.000158+1.0000001 
0.000158  +  1.0000001 
0.000158  +  1.0000001 
0.000158  +  0.9999941 
0.000158  +  0.9999801 
0.000158  +  0.9999941 
0.051051  +  0.1231061 
0.051031  +  0.1231551 

0.1352917  +  0.1151021 
0.134745  +  0.1147821 
0.134744  +  0.1147821 
0.127835  +  0.1234731 
0.127839  +  0.1238441 
0.127646  +  0.1236801 
0.391769  +  0.1604091 
0.390212  +  0.1610581 

0.015698  +  1.000601 i 
0.015702  +  1. 000595 i 
0.015702  +  1.000595i 
0.015784  +  1.000040i 
0.015855  +  0.9985771 
0.015787  +  0.9999911 
0.162270  +  0.387268i 
0.161632  +  0.3888261 

3.863535  +  1.6251971 
3.863438  +  1 . 62541 9i 
3.863425  +  1.6254141 
3.847469  +  1.6329471 
3.862934  +  1.6267131 
3.847456  +  1.6329781 
3.917688  +  1. 604092 i 
3.902121  +  1 .161058i 

1.598486  +  3.928208i 
1.598720  +  3.9281151 
1.598715  +  3.9281031 
1.593312  +  3.9431781 
1.599938  +  3.9275801 
1.593343  +  3.9431651 
1.622698  +  3.8727681 
1.616319  +  3.8882621 
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Table  8  Characteristic  Roots  (w  =  e^cosng) 
Graphlte/Epox 


Equation 


P1’P2 


1.986613 
1 .987408 
1.987389 
1.984245 
1.987400 
1  .987387 
1.987400 
1.987387 


2.111661 

2.119304 

2.119284 

2.024353 

2.120223 

2.027557 

2.120533 

2.027929 


19.925442 

19.925560 

19.925434 

19.889481 

19.925642 
19.889795 

19.925643 
19.889797 


+  1 . 9881861 
+  1 .9889761 
+  1.9889701 
+  1.9905351 
+  1.987400i 
+  1.9873871 
+  1.9874001 
+  1.9873871 


+  1.8548321 
+  1.8635291 
+  1 .8635231* 
+  1.9510961 
+  1.8628971 
+  1.9480141 
+  1.8633581 
+  1.9484341 


+  1 9 .8224531 
+  19.822572i 
+  1 9 . 8224481 
+  19.8582861 
+  1 9.8224961 
+  19.857972i 
+  19.8224981 
+  19.8579731 


0.042075  +  0.036972i 


0.01 3807 
0.01 3807 
0.013807 
0.013791 
0.01 3815 
0.01 3791 
0.015953 
0.015924 


+  0.0137531 
+  0.01 37531 
+  0.0137531 
+  0.0137691 
+  0.0137441 
+  0.0137691 
+  0.0158701 
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